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ABSTRACT 

The synchrotron shock model (SSM) for gamma-ray burst emission makes 
a testable prediction: that the observed low-energy power-law photon number 
spectral index cannot exceed —2/3 (where the photon model is dehned with 
a positive index: dN/dE oc U"). We have collected time-resolved spectral £t 
parameters for over 100 bright bursts observed by the Burst And Transient 
Source Experiment on board the Compton Gamma Ray Observatory. Using 
this database, we hnd 23 bursts in which the spectral index limit of the SSM 
is violated. We discuss elements of the analysis methodology that affect the 
robustness of this result, as well as some of the escape hatches left for the SSM 
by theory. 


Subject headings: gamma rays: bursts — radiation mechanisms: non-thermal 


1. Introduction 


Gamma-ray bursts (GRBs) have been an astronomical puzzle for 30 years. Details 
of the gamma-ray emission from bursts are still elusive, although there is an emerging 
agreement that persistent afterglow emission at other wavelengths from locations consistent 
with those of several recent GRBs can be related to expanding hreballs at cosmological 
distances ([Goodman 1986|, [Metzger et al. 1997|, [Waxman 1997[). One promising mechanism 
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that has been proposed for gamma-ray emission is the synchrotron shock model (SSM): 
synchrotron emission from particles accelerated in a relativistic shock is Lorentz-boosted 
into the gamma-ray band ( |Kees fc Meszaros 1992| , [IVfeszaros fc Rees 1993| , [Rees &: Meszaros 
1994, Katz 191?^, havani 1996). The SSM identihes the spectral break observed in most 


burst spectra with the characteristic synchrotron energy in the emitter’s rest frame, boosted 
into the observer’s frame. Thus, the htted energy of the spectral break would contain 
information about both the bulk Lorentz motion of the emitters and the equipartition 
magnetic held producing the synchrotron radiation (|Tavani 1995|) . In addition, the SSM 
makes the specihc prediction that the low-energy power-law spectral index of the observed 
photon number spectra should not exceed —2/3, with the assumption that the optical depth 
of the shocked material is less than unity. We use the convention that the power-law index 
a takes a positive sign: dN/dE oc E°‘. The —2/3 value is derived from the synchrotron 
single-particle emission spectrum: bulk relativistic motion ensures that the mean particle 
energy is large enough that there are few low-energy particles ( |Katz 19944 [havani 1995|) . 
Below the cyclotron fundamental energy, the low-energy spectrum is approximately that 
produced by mono-energetic particles and is thus independent of their actual distribution 
and should be a constant —2/3 spectral index power law. If that were the total story, 
the model would already be rejected, since it is well-known that bursts are observed 
to have a variety of spectral behaviors at low energies ([Band et ah 199"^ [Preece et ah 


19961 , jStrohmayer et al. 1997] , [Crider et ah 1997[ ). However, it also has been noted that 


the timescale for synchrotron cooling of the particles may be shorter than the duration 
of burst pulses ( [Katz 1994b[ , [Sari, INarayan fc Piran 19^ , [Sari fc Piran 1997[) . A cooling 
distribution of particles is characterized by a power-law index of —2, which translates into a 
—3/2 photon number index through the synchrotron power-law emission formula ([Rybicki 


[fc Lightman 1979[) . If one includes the effects of cooling of the particle distribution, the 
low-energy spectral index can encompass the range of —3/2 to —2/3. However, the spectral 
slope still cannot be greater than the fundamental single-particle limit of —2/3. 

In this Letter, we test the SSM limit on spectral behavior by examining how well the 
data support it. We draw our results from a catalog of time-sequences of spectral hts to 137 
bursts selected for their high flux and fluence, using (mostly) Large-Area Detector (LAD) 
data from the Burst And Transient Source Experiment (BATSE) on the Compton Gamma 
Ray Observatory (see Preece et al. 199'8| for details of the analysis). In the next section, we 
determine the best spectral indicators derivable from the data and show how they should 
bracket the ‘true’ low-energy spectral behavior. In §3, for each of the bursts in the catalog 
we compare the effective low-energy power-law spectral index with the SSM limit line. 
Finally, in §4 we discuss the implications these results have on further theoretical modeling. 
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2. Spectral Modeling 


The model typically chosen for spectral htting of bursts is the empirical ‘GRB’ function 
( IBand et al. 1^^ ): 
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where the four model parameters are: the amplitude A, a low-energy spectral index a, 
a high-energy spectral index f3 and an energy -Epeak that corresponds to the peak of the 
spectrum in vTy if (3 is less than —2. In this expression, -Epeak and a are jointly determined 
by the low-energy continuum, while (3 is solely determined by the spectrum above the break 
energy, E^break- If this break energy lies above the highest energy available to the detector, 
(3 is ill-dehned, so we must substitute a related form of the model with the high-energy 
power-law omitted. In cases where a sharp inherent curvature of the spectrum results in 
an unacceptable value of for a £t to the GRB spectral form, a broken power-law (BPL) 
model often generates better results. 


Observed burst spectra are such that the data rarely approach the GRB spectral 
low-energy power-law within the energy range of BATSE and most other burst experiments. 
At the low end of the spectrum, a is approached only asymptotically, as seen by comparing 
the dotted line in Figure |l] with the GRB model £t to an example spectrum {solid line). 
The model-dependent photon ‘data’ rates consist of photon model rates weighted by the 
ratios of the deconvolved model rates and the count rates in each data channel, allowing 
different photon models to be represented simply on a single plot. A better measure of the 
actual lower-energy behavior is an effective spectral index, e. g., the slope of the power-law 
tangent to the GRB function at some chosen energy (Efid), which can be found analytically: 
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For typical values of a, i. e., a > —2, Q;eff(Efid) will always be less than a by an amount 
that depends upon the value chosen for Efid- We chose Em = 25 keV for the following 
reasons: The observed spectrum is consistent with the GRB function above about 25 keV; 
below this value, deviations from the standard GRB function have been observed that may 
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indicate the existence of an additional low-energy spectral component ( Preece et al. 19^ ). 
In addition, 25 keV is just greater than the low-energy cut-off for the LADs; it is above 
an electronic spectral distortion in the SDs at typical gain settings ( [Band et al. 19921 ) and 
it is also below htted values for -Epeak in nearly all bursts. The effective spectral index at 
25 keV will be denoted Q! 25 - Once has been hxed, the correction gets smaller as the 
htted value for T^peak gets larger, since is closer to the asymptotic value a. This is also 
seen in Figure |I], where F^peak = 1308 keV, so a = —0.04 Ri Q !25 (the GRB function is not 
appropriate for this spectrum and is shown for illustration only). In any practical case, 025 
serves as an upper bound on the low-energy spectral index, in the sense that it is more 
positive than an average slope through the model. We will use 0:25 in this work, rather than 
a, since it is more indicative of the slope actually seen in the data. 


To the extent that all burst spectra have continuous curvature rather than a sharp 
break, a BPL function £t (Figure - dashed line) will generate a spectral index (which we 
will call Q!pl) that is a lower bound on the ‘true’ low-energy spectral behavior, in that it will 
be more negative than it would be for a model that takes into account some curvature, such 
as the GRB function. Thus, to be conservative in our estimate for how well the BATSE 
data support the SSM, app provides the best available comparison, since if it lies above the 
limiting —2/3 line, we can be reasonably sure the ‘true’ spectrum will violate the SSM. Of 
course, this will do us no good if the BPL model is not an acceptable fit to the data. Thus, 
we use hts to two models to put bounds on the ‘true’ low-energy spectral behavior: a is a 
model parameter asymptotic to the data, but we can approximate the upper bound with 
the quantity 0:25 dehned above. The BPL estimate, app, is a lower bound. 


3. Observations 

In Figure we show results from our BATSE catalog of time-resolved spectroscopy of 
bright bursts (peak hux > 10 ph cm“^ s“^ and huence > 4 x 10“^ erg cm“^), using mostly 
LAD data. The value of 025 is plotted against Epeak for the spectrum in each burst for 
which 025 reaches its maximum value. In addition to the 0:25 ^ Epeak parameters obtained 
by fitting the GRB function, we have also included opp - Epreak pairs (diamonds) for those 
bursts where use of the BPL model was more appropriate. Overlaying this plot is the 
synchrotron —2/3 spectral index ‘death line’ [dashes), as well as the lower bound allowed for 
cooling distributions [dotted line). The SSM is quite safe within these boundaries; however, 
44% of the total have maximum low-energy spectral indices above —2/3, in a region that 
is totally excluded in the SSM. In particular, 32% of the BPL bursts have points that he 
above this line. With some correlation between the two displayed parameters, the most 
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stringent limits are set by the spectral indices of spectra with the highest values of -Epeak, 
as well as any bursts that have a BPL index greater than —2/3. Clearly, many such bursts 
exist; one is the example shown in Figure |^, which is a BPL £t with apL = —0.263 and 
Ehreak = 456 keV, the highest-energy diamond above the —2/3 limiting line. 

The error bars shown in Figure ^ represent the la error for each parameter considered 
singly, as obtained from the covariance matrix of the £t. If the errors in the determination 
of the spectral index were not normally distributed, perhaps having a broad tail instead, 
then values that are much greater than the SSM limit of —2/3 might violate the SSM. 

To test this, we created 4 sets each of 1000 simulated spectra, based on the GRB model 
spectrum propagated through the detector response, with random Poisson fluctuations 
determined for the counts in each data channel. These consisted of bright, medium and dim 
spectra, with expected errors of 0.06, 0.2 and 0.34, respectively, for the htted value of a. 
The parameter values assumed were Ebreak ~ 400 keV, (5 = —2 and a = —2/3, to provide 
a worst case test. A second set of medium-brightness spectra was also created, to test the 
effect of a small assumed value of Ebreak = 200 keV. A histogram of the resulting htted 
a values is approximately consistent with a Gaussian distribution for each set, however, 
distributions for the dimmer sets of spectra have slightly extended tails on the side of 
more positive values of a. These larger values all have large computed errors associated 
with them, so the resulting distributions of deviations of a from the mean, weighted by 
their errors (see below), do not have a tail. Rather, the average over the 4 simulations 
for the number of error-weighted deviations greater than 2 a is 7.25, compared with the 
22.7 deviations expected for 1000 trials in a normal distribution. Interestingly, the effect 
where large values of a have larger associated errors can be observed in Figure where the 
errors above the —2/3 line are, on average, 1.5 times those below. The KS probability we 
calculate from comparing the a > —2/3 tails between our observations and the simulations 
is 3 X 10“^® (Dks = 0.26 for A^^obs = 1012 and = 1689), indicating that it is quite 
unlikely that the two distributions are the same. 

Another way to view these results is as a histogram of = {a — (—2/3))/(Tq (here, 
we use a to stand in as a generic low-energy power-law index); that is, the deviation of a 
from —2/3 in units of the standard error for each £t. Figure ^ shows the distribution of Aa^ 
for the entire ensemble of 3957 spectra from the 137 bursts £t, where the SSM-violating 
region consists of the positive values to the right of zero. There are 312 spectra total in the 
bins containing 3cr and greater (that is: > 2.4cr). This may be compared with our expected 
error distribution, based upon the 4000 total simulated spectral hts described above (Fig. ^ 
- dashed line), where the total over the same bins is 8 (in fact, 0 > 3.4cr). We may estimate 
the probability that 312 out of 3957 spectra would arise from 2.4 ct huctuations larger than 
—2/3, assuming the SSM prediction that no such spectra exist. The simulations show that 
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the chance probability of a —2/3 slope spectra having a £t value more than +2.4(T away 
from —2/3 is ~ 8/4000. Using the binomial probability distribution, the chance that the 
observed 312 out of 3957 could result from random fluctuations is < 10“^®®, showing the 
observations to be quite inconsistent with the SSM. 


4. Discussion 

We have shown that there are a large number of bursts that violate the limit on 
low-energy spectral behavior imposed by the basic synchrotron emission mechanism acting 
in a relativistic shock. It is worthwhile considering which spectral models can accommodate 
this observation. First of all, it has been suggested ([Liang et ah 19971) that Compton 
upscattering of soft photons by an energetic distribution of particles can signihcantly 
modify the basic synchrotron emission spectrum, with energetic particles boosting their own 
synchrotron emission into the observed gamma-ray band. The details of the spectral shape 
thus depend upon the particle distribution, as well as the shape of the photon spectrum 
at the energies where it is being sampled for upscattering. The observed steep low-energy 
spectral indices would arise for certain combinations of the source parameters and then 
would evolve to smaller values as the particles cool. An important prediction from this 
model is that a low-energy component, independent of the observed gamma-ray emission, 
must be present with sufficient strength to serve as a pool of photons for upscattering. It is 
indeed possible that this separate component may already have been observed (jPreece et al. 

. If the low-energy portion is truly an independent component, 
it will have an independent time history as well. The troubling part of this idea is that all 
bursts should have this component to some extent, since there is no evident bi-modality to 
the low-energy behavior (cf. Figure |^) that would indicate that the SSM-violating bursts 
are somehow different. 

Synchrotron self-absorption (SSA) is another mechanism that would tend to increase 
the low-energy continuum spectral index. The maximum photon spectral index that could 
be observed is -1-3/2, so all of the bursts presented here would be consistent. The photon 
opacity must be close to unity in all these sources for some to be self-absorbed. For SSA 
to work the optical depth must be greater than one at energies below Upeak- Thus, we 
trade one mystery for another: the narrow distribution of Upeak arises in the fact that many 
bursts have optical depths close to unity, rather than from a narrow distribution of Lorentz 
factors (see next paragraph). Also, if the photon density is high, it may be very difficult to 
overcome the opacity arising from photon-photon pair-production. 

Since the observed value of Upeak should scale as the bulk Lorentz factor of the emitting 


1996, Strohmayer et ah 1997 
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material to the fourth power, a narrow observed distribution of -Epeak implies either that the 
rest-frame value is extraordinarily precisely determined for all bursts or that the Lorentz 
factors of the entire ensemble lie in a very narrow distribution. Since neither alternative 
is satisfactory, Brainerd (1994)| proposed that the observed spectra may arise from much 
different spectra that have been attenuated by intervening material through relativistic 
Compton scattering. The low-energy behavior is directly related to conditions at the source, 
namely the optical depth of the scatterers. Large optical depths will result in steep htted 
low-energy power-law indices, and the distribution of these should somehow relate to the 
distribution of densities of the material surrounding putative sources. 


Thanks to B. Schaefer who suggested adding error bars to the a - Ep^ak figure from 
Preece et ah 1996. The anonymous referee has contributed suggestions that have lead to 
an improved error analysis. 
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Fig. 1.— Photon spectrum accumulated between 0.448 - 0.768 s after the trigger for 
3B910814. Two best-£t models are plotted on top of the deconvolved data: the GRB function 
{solid line) and the BPL model {dashed line; apL = —0.263). The tangent slope at 25 keV 
is also shown: a = —0.04 {dotted line). Neither model is as steep as —2/3 at low energies. 
The errors represent model variances. 
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Fig. 2.— A plot of fitted low-energy power-law index against E^^ak- The spectrum chosen 
for each burst is the one with the largest value of a. Diamonds indicate bursts £t with a 
broken power-law model. The SSM-violating region is bounded from below by the ‘death 
line’ {dashed) and the acceptable region is bounded by the —3/2 line {dotted). Points where 
the error bars on 0:25 exceed the plot area were omitted for clarity. 
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Fig. 3.— A histogram of deviations of the low-energy power-law indices from —2/3 in units 
of their la error for 3957 htted spectra. Positive deviations represent a violation of the SSM. 
Also plotted is the distribution of la deviations from the mean for spectra indices obtained 
from hts to 4000 simulated spectra {dotted). The peak has been clipped so that the wings 
of the two distributions can be compared. 







